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The kinetics of chiral transitions in quark matter is studied in a two flavor Nambu- 
Jona-Lasinio model. We focus on the phase ordering dynamics subsequent to a 

o , 

^SJ , temperature quench from the massless quark phase to the massive quark phase. We 

' study the dynamics by considering a phenomenological model (Ginzburg-Landau 
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free-energy functional). The morphology of the ordering system is characterized by 
the scaling of the order-parameter correlation function. 



Heavy-ion collision experiments at high energies produce hot and dense strong! 
D ■ interacting matter and provide the opportunity to explore the phase diagram of QCD 



It is worthwhile to stress that, in a phase transition process, information about which equi- 



^ \ librium phase has lowest free energy is not sufficient to discuss all possible structures that 
Tjj- ■ the system can have. One has to understand the kinetics of the process by which the phase 

00 ; 

, ordering or disordering proceeds and the nature of nonequilibrium structures that the sys- 

tem must go through on its way to reach equilibrium J^, ^. In the present work, we study 



kinetics of chiral transitions in hot and dense matter 



We use the two-flavor Nambu-Jona-Lasinio (NJL) model to study the chiral symmetry 
^ ■ breaking in QCD j^. The expression for the thermodynamic potential is obtained in the 



mean-field approximation: 



n(M, /3, ^^) = --^ fdk \\n\l + e-/^(v^™"-M)l + In [l + ^-P{VWTW+,)] | 
{2n)-^p J ^ I J L J 

dk (\/k^ + M^ -\k\) + (1) 
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Here, we have taken vanishing current quark mass, and introduce the constituent mass 
M = —2gps with = {ijjip) being the scalar density and g = G[l + l/{4:Nc)]. The details 
of the mean-field approximation are reported elsewhere Q]. 
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The phase diagram resulting from Eq. ([T]) is shown in Fig. [T]^a). For the parameters of 
the NJL model, we have taken the three-momentum cutoff A = 653.30 MeV, and coupling 
G = 5.0163 X 10~^ MeV~^ With these parameters, the vacuum mass of quarks is 
M ~ 312 MeV. At T = 0, a first-order transition takes place at ;U ~ 326.321 MeV. For 
/i = 0, a second-order transition takes place at T ~ 190 MeV. The first-order line (I) meets 
the second-order line (II) at the tricritical point {f^tcpjTtcp) — (282.58,78) MeV. 

Close to the phase boundary, the potential in Eq. ([1]) may be expanded as a Ginzburg- 
Landau (GL) potential in the order parameter M: 

n{M) = n{0) + -M'' + -M^ + -M^ + OiM''), (2) 
2 4 6 



correct upto logarithmic corrections [8|. In the following, we consider the expansion of 
potential Q (M) upto the M^-term. The first two coefficients in Eq. ([2]) can then be exactly 
obtained by comparison with Eq. ([T]) while the coefficients b and d are obtained by fitting 
the GL potential with the potential Eq. ([T]) for a given T and /x. For stability, we require 
d>0. 

The extrema of the potential in Eq. ([2]) are determined by the gap equation Cl'{M) = 
with solution M = 0, and M| = {-b ± y/W^^Aad) / {2d) . For b > 0, the transition is 
second-order, with the stationary points M = (for a > 0) or M = 0, ±M_|_ (for a < 0). 
For a < 0, the preferred equilibrium state is the one with massive quarks. For 6 < 0, the 
solutions of the gap equation are as follows: (i) M = for a > b"^ /{Ad), (ii) M = 0, ±M+, 
±M_ for 6^/(4(i) > a > 0, and (iii) M = 0, ±M+ for a < 0. A first-order transition takes 
place at = 3b'^/{16d) with the order parameter jumping discontinuously from M = to 
M = ±M+ = ±(3|6|/4d)i/2. 

The evolution of the system is described by the time-dependent Ginzburg-Landau 
(TDGL) equation: 

which models the over-damped relaxational dynamics of M(r, t) to the minimum of Q [M] 
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Here, F is the inverse damping coefficient, and 9{f,t) is the noise term satisfy- 
ing the fluctuation-dissipation relation: {9{f,t)) = and (9{r' ,t')9{r" ,1")"^ = 2TT6{r' — 
r")5{t' — t"). We use the natural scales of order parameter, space and time to introduce 



dimensionless variables: M = MqM' (Mq = ^\a\/\b\)] f = ^r' (^ = ^JK/\a\)■, t = rt' 
[t = (r|a|)^^]; 9 = (r|a|^/^T^/^/|6|^/^) 9'. Dropping the primes, we obtain the dimensionless 
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TDGL equation: 
(9 

— M (f, t) = -sgn (a) M - sgn (6) - XM^ + V^M + 9 (f, t) , (4) 

L/(/ 

where A = |a|rf/6^ > 0. 

First, we study the ordering dynamics for 6 > corresponding to a quenching through the 
second order transition. We solve Eq. (jlj) numerically using an Euler-discretization scheme 
with an isotropic Laplacian implemented on a d=2 lattice of size A^^ (A^ = 4096), with 
periodic boundary conditions in all directions. The details of the simulation can be found 
elsewhere {g]. In Fig. |2l we show the evolution of a disordered initial condition for Eq. (jl]) 
with 6 > and a < - corresponding to a temperature quench through the second-order 
line (II) in Fig.[T]^a). The initial state consists of small-amplitude thermal fluctuations about 
the massless phase M = 0. The system rapidly evolves into domains of the massive phase 
with M ~ M+ and M ~ — Af_|_. The interfaces between these domains correspond to M = 
- their evolution is shown in the snapshots (frames at the top) of Fig. |2l The frames at the 
bottom show the order parameter variation of the snapshots at the top. 

The domains have a characteristic length scale L{t), which grows with time. The growth 
process is analogous to coarsening in the TDGL equation with an M^-potential {2!, yj. The 
order-parameter correlation function C{r,t) shows dynamical scaling C{r,t) = f{r/L). The 
scaling function f{x) = (2/7r) sin^^(e~^^) has been calculated by Ohta et al. (OJK) in 
the context of an ordering ferromagnet. Further, the domain scale obeys the Allen-Cahn 
(AC) growth law, L{t) ~ t^s xhe same growth law has also been obtained via a closed 

;ime path formalism of relativistic finite-temperature field theory applied to the NJL model 

III- 

Next, we consider the case with b < . In this case, a first-order chiral transition occurs 
for a < ac = 3b'^/{16d) (or A < Ac = 3/16). For a < 0, the potential has a double-well 
structure and the ordering dynamics is equivalent to M*^-theory, i.e., the domain growth 
scenario is similar to Fig. |2l We focus on a quench from the disordered state (with M = 0) 
to < a < or < A < Ac, corresponding to a quench between the first-order line (I) and 
Si in Fig. [TJ^a) - the corresponding points are denoted by crosses. The massless state is now 
a metastable state of the M^-potential. The chiral transition proceeds via the nucleation 
and growth of droplets of the preferred phase (M = ±M+). The nucleation results from 
large fiuctuations in the initial condition or thermal fluctuations during the evolution. In 
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Fig. [3l we show the nucleation and growth process. At early times (t = 400), the system 
is primarily in the M = phase with small droplets of the preferred phase. These droplets 
grow in time and coalesce into domains. 

In the late stages of growth, there is no memory of the nucleation dynamics which char- 
acterized growth during the early stages. Let us note that a large dissipation (i.e., F"^) will 
make the equilibration time larger. If this becomes larger than the life time of the fire ball 
in heavy ion collision, the system may linger in the symmetry restored phase longer even 
when the temperature has already decreased below Tg. This means the calculated based 
on equilibrium thermodynamic models could be higher than the value that shows up in the 
experiments 0, 6|. 

To summarize: we have studied the kinetics of chiral phase transitions in QCD. In terms 
of the quark degrees of freedom, the phase diagram is obtained using the NJL model. An 
equivalent coarse-grained description is obtained from an M^-Ginzburg-Landau (GL) free 
energy. The chiral kinetics is modeled via the TDGL equation, and we consider the over- 
damped case. We study the ordering dynamics resulting from a sudden temperature quench 
through the first-order (I) or second-order (II) transition lines in Figs. HJ^a). For shallow 
quenches through II and deep quenches through I, the massless phase is spontaneously un- 
stable and evolves to the massive phase via spinodal decomposition. For shallow quenches 
through I, the massless phase is metastable and the phase transition proceeds via the nucle- 
ation and growth of droplets of the massive phase. The merger of these droplets results in 
late-stage domain growth analogous to that for the unstable case. In all cases, the asymp- 
totic growth process exhibits dynamical scaling, and the growth law is L(t) ~ t^/^. Given the 
dynamical universality of the processes involved, our results are of much wider applicability 
than the simple NJL Hamiltonian considered here. 
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FIG. 1: (a) Phase diagram of the Nambu-Jona-Lasmio (NJL) model m the (/i,T)-plane for zero current 
quark mass. A hue of first-order transitions (I, green online) meets a line of second-order transitions (II, 
blue online) at the tricritical point (tcp). We have [ntcp^Ttcp] — (282.58, 78) MeV. The dot-dashed lines 5*1 
and 5*2 denote the spinodals or metastability limits for the first-order transitions. The open symbols denote 
4 combinations of (/i,T), chosen to represent qualitatively different shapes of the NJL potential. The cross 
denotes the point at which we quench the system for 6 < 0. (b) Plot of f2 (Af, /3, /i) from Eq. ([T]) as a function 
of M. The (/i, T)- values are marked in (a). The solid lines superposed on the potentials correspond to the 
GL potential in Eq. ([2]). 
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FIG. 2: Interface evolution after temperature quench through second-order hue (II) in Figs. [ija). The 
snapshots at the top show the interfaces (A/ = 0) at t = 200, 4000 (in units of r). They were obtained by 
numericahy solving Eq. Q as described in the text with a < 0, 6 > 0, A = 0.14. The noise amplitude was 
e = 0.008. The frames at the bottom show the variation of the order parameter along the diagonal. 
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FIG. 3: Interface evolution after temperature quench through first-order line (I) in Figs.[lja). The snapshots 
at the top show the interfaces (A/ = 0) at t = 400, 4000. They were obtained by solving Eq. Q with 6 < 0, 
Cc > a > and A = 0.14. The frames at the bottom show the variation of the order parameter along the 
diagonal. Notice that the metastable patches (M — 0) at i = 400 are absent at later times. 



